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Abstract: We study extremal black hole solutions to four dimensional N = 2 supergra- 
vity based on a cubic symmetric scalar manifold. Using the coset construction available 
for these models, we define the first order flow equations implied by the corresponding 
nilpotency conditions on the three-dimensional scalar momenta for the composite non-BPS 
class of multi-centre black holes. As an application, we directly solve these equations for 
the single-centre subclass, and write the general solution in a manifestly duality covariant 
form. This includes all single-centre under-rotating non-BPS solutions, as well as their 
non- interacting multi-centre generalisations. 
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1. Introduction and Overview 

Black holes in the context of string theory have been a long lasting field of research, 
due to its deep connection to fundamental aspects of the theory. One important facet of 
the subject has been the construction of supergravity solutions describing the low energy 
strong coupling regime of these systems. For super symmetric black holes, it has been 
possible to find these solutions explicitly, making use of the constraint imposed by the 
residual supercharges, even including higher derivative corrections, see e.g. (l|, |2], [| [|, || 
H [7|, H [9|. Such a task however has proven more difficult for more general black holes, 
since one is forced to consider the full second order equations of motion rather than the 
first order BPS conditions. The simplest generalisation is the class of extremal black holes, 
which are still characterised by a vanishing Hawking temperature, but do not preserve 
any super symmetry. The corresponding static solutions are known to be described by first 



order equations as well M, O, P, pL 14, 15, 16, 17, IS, IS, GO, 21], although the latter 



are then not a direct consequence of supersymmetry. 

Extremal black holes in supergravity theories fall in two distinct categories, namely 
the under-rotating and the over-rotating branches, where the former contains several sub- 
classes. The over-rotating (or ergo) branch is characterised by the presence of an ergo-region 
and includes the extremal Kerr solution. In contrast, we will focus on the under-rotating 
(or ergo-free) black holes, which then admit a flat three-dimensional base, and include the 
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24]. Single centre under-rotating non-BPS black holes 



static extremal black holes [p2j , 
have been studied throughout the last decade or so, from various aspects and using various 
techniques, see for example [25, 26, 27, 28, 10, 15, pi], 12, 17, 15] and references therein 
for some developments. Using the seed solution of 11, 13, 29] combined with a general 
duality transformation as explained in [30], one can construct any desired solution, but a 
manifestly duality covariant formulation was lacking. 

For theories coupled to a symmetric scalar manifold, as the ones we will deal with, 
there are three classes of solutions describing interacting ergo- free extremal black holes, 
distinguished by the algebraic properties of the first order systems that describe them. 
These are the standard BPS solutions [0, ||, the so-called almost BPS solutions [31], and 
the composite non-BPS solutions ||3"2}] . After the BPS class, the almost-BPS class is perhaps 
the best studied, especially in the five dimensional uplift, where it was originally discovered 
and extended, see e.g. pil p3, El 



p3| , [34| , pq , pq [. On the other hand, the composite non-BPS class 
appears to be simpler in the four dimensional context and is the class we discuss in the 
following. 

In this paper we will consider the equations of motion for stationary solutions as 
described by the dimensional reduction of the theory along time, i.e. as a non-linear sigma 
model over a pseudo-Riemannian symmetric space coupled to three-dimensional Euclidean 
gravity J37j. Spherically symmetric black hole solutions correspond to geodesies on this 
symmetric space and, in particular, extremal static black holes correspond to the subclass of 
null geodesies. It was shown in |38|, |l4], ^9|, [5], 41] that null geodesies associated to extremal 
black holes are classified by the nilpotent orbits in which the corresponding Noether charges 
lie in, thus recasting the equations of motion to an eigenvalue equation for the momentum 
of the coset scalars. The classification of extremal solutions in terms of nilpotent orbits 



has been studied in details in [17, 2C, 42 



Considering a general stationary Ansatz with a flat three dimensional base, one finds 
that the only regular solutions are described by scalar fields taking values in a nilpotent 
subgroup of the three-dimensional duality group. Using the nilpotent orbits classification, 
one can determine pertinent nilpotent subalgebras, which lead to the description of station- 
ary solutions describing interacting black holes. It was shown that all the three systems, 
BPS, almost BPS and composite non-BPS can be described in terms of associated nilpotent 
subalgebras [32, 44]. In these papers, explicit examples of these solutions were constructed 
in terms of specific nilpotent subalgebras. 

However, these seed solutions are not duality covariant by construction, whereas one 
should be able to construct them without referring to a specific representative. Although 
the general solutions can be obtained from these seed solutions by duality transformations, 
the resulting form of the solutions is not easy to parametrize in terms of physically rele- 
vant quantities. A duality covariant formulation of these solvable systems would permit 
to obtain these solutions in a form exhibiting their physical properties. This is of particu- 
lar importance when addressing issues such as existence and stability of composite bound 
states. For the BPS solutions (as well as the non-BPS = 0), the first order system asso- 
ciated to the corresponding nilpotent orbits was written in a manifestly duality invariant 
form [45 1, leading to a generalisation of the BPS solutions j|, [| to M = 8 supergravity, 
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and non-BPS = solutions. 

In what follows, we explain how the composite non-BPS solutions can also be described 
in a manifestly duality invariant form. As it turns out, the system is characterised by a 
constant very small charge vector (i.e. such that its quartic invariant satisfies I4 = dl^ = 
c^Zilad = 0), which has only one charge component in an appropriate duality frame. This 
represents an auxiliary variable that restricts the types of charges allowed in the various 
centres, or conversely is fixed by the physical charges for a regular solution. 

As a first application, we will solve explicitly the system in the restricted case describing 
single centre (or non-interacting multi-centre) black holes in a manifestly duality covariant 
way. The reduction to this subclass is effected by a suitable reality constraint on the scalar 
momenta, which we use systematically to simplify the problem. The result is a stabilisation 
equation for the scalars throughout the black hole background that parallels the one derived 
in [Q, U for supersymmetric solutions. However, we find new terms, proportional to the 
very small vector driving the flow and its magnetic dual, with coefficients depending on the 
harmonic functions carrying the electromagnetic charges and angular momentum. These 
represent a duality covariant realisation of the simple change of sign for a particular charge, 
which has been observed to relate some simple non-BPS solutions to supersymmetric ones 
1 26, 27, 46], and the non-harmonic term obtained in [[HJ] by considering the seed solution 



of []29|, respectively. 

Understanding how to solve the first order system when restricted to single centre 
solutions is a necessary step towards the resolution of the composite non-BPS system. 
In order to do so, one must rewrite the non-linear first order system we describe in this 
paper, in a second order linear system of differential equations. The explicit solution of the 
composite non-BPS system will be the purpose of forthcoming research. 

We start by setting up notation and giving some background on both four dimensional 
N = 2 supergravity and the corresponding non-linear sigma model obtained by time-like 
reduction in section |2[ In section || we proceed to the discussion of the algebraic structure of 
the composite non-BPS system and discuss the reality constraint that reduces it to describe 
single centre solutions. We then go on to solve explicitly the flow equations for the single 
centre class and discuss some aspects of the solutions obtained in section |j. We conclude 
in section I with some general remarks and future plans of extending to the multi centre 
classes, whereas the two appendices are devoted to the illustrating the general solution 
with the seed solution of [^9| and to the derivation of the duality invariant constraints on 
the charges. 



2. Non-linear sigma model formulation of stationary solutions 

In this section, we collect various formulae and conventions which will be essential for the 
connection of the objects appearing in the three-dimensional non-linear sigma model formu- 
lation of N = 2 supergravity describing stationary solutions to standard four-dimensional 
supergravity variables. 
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2.1 M = 2 supergravity and symmetric special Kahler geometry 

The bosonic Lagrangian of J\f = 2 supergravity coupled to n v vector multiplets reads 



Svre- 1 /; = -^R-i(D^V,D^V)-lF^G^ , (2.1) 

where the F* v = d^A^ - d u Aj 1 for I = 0, . . . encompass the graviphoton and the gauge 
fields of the vector multiplets and Gj U are the dual field strengths, defined in terms of 
the Fp V though the scalar dependent couplings, whose explicit form will not be relevant in 
what follows. The gauge field equations of motion and Bianchi identities can then be cast 
as a Bianchi identity on the symplectic vector 

^ - (<£) ■ (2 ' 2) 

whose integral over any two-cycle defines the associated electromagnetic charges through 



1 

2^ 



T. (2.3) 

s 2 



The physical scalar fields f, which parametrize a special Kahler space Ai^ of complex 
dimension n v , only appear in (|2.1|) through the section, V, of a holomorphic U(l) x Sp(2n v + 
2, R) bundle over M^. Choosing a basis, this section can be written in components in terms 
of scalars X 1 as 

OF 




>' - ( , J • F i = dxi> ( 2 - 4 ) 

where F is a holomorphic function of degree two, called the prepotential, which we will 
always consider to be cubic 

1 X*XJX k ^ M[X\ 
F = ~Q C ijk x0 =-^pT> ( 2 - 5 ) 

for completely symmetric Cj,^, %= 1, . . . n v , and we introduced the cubic norm jV[X]. The 
section V is subject to the constraint 

(V,V)=i, (2.6) 

and is uniquely determined by the physical scalar fields t % up to a local f/(l) transformation. 
The U(l) gauge invariance of fl2.l| ) is ensured by the appearance of the Kahler connection 
in the covariant derivative. The Kahler potential on .M4 is defined up to an arbitrary 
holomorphic function f(t) as 

JC = -ln(iN[t-t\)+f(t) + f(i) (2.7) 

and we fixed the U{\) gauge invariance in terms of Kahler transformations by requiring 
that the Kahler connection is determined by the Kahler potential as 

Q = lm [diKdt 1 ] , (2.8) 



- 4 - 



such that 

g i3 = didjK , D^V = {d^ + i Q M )V = AV = {d t V + \dJC V) dj , (2.9) 

where D{V is the corresponding Kahler covariant derivative on the components of the 
section. With the prepotential (2.5), the special geometry identities [49] reduce to 

DjDiV = g i3 V , DiDjV = ie K c ijk g k ~ k D- k V , (2.10) 

which will be used extensively in what follows. 

We introduce the following notation for any symplectic vector J 

Z(J) = (J,V), (2.11) 
Z i (J) = {J,D i V), (2.12) 

with the understanding that when the argument is form valued, the operation is applied 
component wise. For instance, the central charge of the gauge field is 

Z(F) = ef (G + ?Gi + \ci jk MF k - M[t]F°) , (2.13) 

for the prepotential ( |2.5|) . With these definitions it is possible to introduce a scalar depen- 
dent complex basis for symplectic vectors, given by (V, -DjV), so that any vector J can be 
expanded as 

J = 21m [-Z(J) V + g lj D,Z{J) DjV) , (2.14) 
whereas the symplectic inner product can be expressed as 

(J U J 2 ) = 21m [-Z(Ji) Z{J 2 ) + Z a (J t ) Z a (J 2 )} . (2.15) 



Finally, we introduce the notion of complex self duality of the gauge fields ( |2.2|) , which 
satisfy the identity 

JJ- = _*j- 5 (2.16) 

where J is a scalar dependent complex structure defined as 

JV = -iV , JAV = iDiV ■ (2.17) 

In this paper we will consider that M4 is moreover a symmetric space, such that the 
coefficients are left invariant by the action of a group G5. In this case one can define 
the vielbeins 011M4 such that 

9 iJ = eie a \ (2.18) 

and 

C-abc = i S. C a C^C c Cjjfc) (2-19) 

where c a b c is a constant symmetric tensor left invariant by the action of K^, which is the 
compact real form of G5. Then, the contravariant symmetric tensor c abc in the conjugate 
representation satisfies the Jordan identity [^] 

Cf(abCcd)gC eI9 = o 6 (a C bcd) ■ (2-20) 
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The symmetric space M.4 is defined as the coset space of the four-dimensional duality 
group G4 by its holonomy subgroup 

M A = (U(l) x K 4 )\G 4 . (2.21) 

The scalar fields can then equivalently be described by a coset representative v in G4, and 
the associated Maurer-Cartan form in the coset component C n " = 04©(u(l)©{4) is defined 
as 

dvv~ l + W = e^df Y a + e aT df l Y a (2.22) 

where Y a (and Hermitian conjugate Y a ) define a basis in g4© (u(l) ©64). Combined with 
( |2,9| ), this equation permits to relate the expression of the Maurer-Cartan form to the 
derivative of the section V in the tangent frame. Similarly, we also define 

Z a (F) = eiDiZ(T) = (F, e l a AV) , (2.23) 

which transforms in the same C n " representation of K4. 

2.2 Time-like reduction and para-quaternionic geometry 

In order to describe stationary asymptotically flat extremal black holes, we introduce the 
standard Ansatz for the metric 

ds 2 = -e 2U (dt + u) 2 + e~ 2U dx ■ dx , (2.24) 

in terms of a scale function U(x) and the Kaluza-Klein one-form oj{x) (with spatial com- 
ponents only), which are both required to asymptote to zero at spatial infinity. Here and 
henceforth, all quantities are independent of time, so that all scalars and forms are defined 
on the flat three-dimensional base. The gauge fields are decomposed in a similar fashion 
as 

2A = C(dt + v) + w (2.25) 
and accordingly for the field strengths 

2F = d( (dt + w) + F , F = ( duj + dw , (2.26) 

where we defined the gauge field scalars arising as the time component of the gauge fields, 
and the one-forms w describing the charges. Here, F is defined as the spatial component 
of the field strength, and is not closed but satisfies 

dF = e 2U * duj A JF , (2.27) 

according to ( 2.16| ), which reads 

d( = e 2U J*F. (2.28) 

Note that this first order equation determines the £ in terms of the vector fields w and the 
scalars. 

The scalar field a dual to the Kaluza-Klein vector to 

e iU *duj = da + {(,d(), (2.29) 
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defines the coordinate of an S 1 fiber over the symplectic torus T parametrized by the 
£'s. Altogether with the scaling factor U and the moduli t l , these fields parametrize the 
para-quaternionic symmetric space 1 

M 3 = G 3 /[SL(2) x Gi) . (2.30) 

This defines the so-called c*-map, which can be related to the standard c-map 51] by 
analytic continuation. The three-dimensional symmetry group Lie algebra 03 decomposes 
as 

03 i(- 2 > e © (&h © 4 ) (o) e © i (2) , (2.31) 

where the weights refer to the eigenvalues under the adjoint action of the 0(i generator. 
The grade one generators in are associated to the gauge invariance with respect to a 
constant shift of the scalars £, and accordingly the grade two generator correspond to the 
invariance under shift of the scalar a. 

At this stage it is important to introduce some properties of the 03 algebra. The 
components of an element of the Lie algebra 03 in the coset component 03©(sl2©04), can be 
decomposed in terms of U(V) x K4 irreducible representations as two complex parameters w 
and Z and two complex vectors Z a , E a which transform in the G nv representation of K4 (the 
same as the scalar fields momenta c^dt 2 in four dimensions) such that A/"[^] = gC abc 
is K4 invariant. The quadratic trace invariant defines the SL{2) x G4 invariant norm 

|w| 2 -|Z| 2 -Z a Z a + £ a £ a , (2.32) 

and the 5X2 algebra is realised on these components as 

5w = ipw + XZ, 5Z = -ipZ + Xw, 5Z a = i P Z a + XT, a , <5£ a = -ipY> a + XZ a , (2.33) 

for a complex A. Closure of this algebra can conveniently be checked in the Cartan complex, 
considering anticommuting parameters A, p and the differential 5 

bp = iXX , 5X = -2ipX , (2.34) 

whose nilpotency is equivalent to the Jacobi identity, similar to the BRST formalism in 
gauge theories. In the same way, the 04 algebra is realised in terms of the elements of £4, 
denoted by G a b, a real 7 and a complex vector A , associated to the decomposition 

04 = u(l) © U e C n " , (2.35) 

which defines the symmetric space A^4. The action of 04 on the coset component can be 
written as 

5w = A a Z a + 3i 7 w 5£ a = A a Z + c abc A b Z c + G a fe S fe + i 7 S a , 

8Z = A a S a + 3ijZ 8Z a = h a w + c abc h b ^ c + G a b Z b + i 1 Z a . '~' ><)! 



1 Here para-quaternionic refers to the property that the holonomy group of M3 SL(2) x G4 C SL(2) x 
Sp(2n v ,R). 
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The corresponding algebra is realised in terms of anticommuting parameters with the nilpo- 
tent differential 

5A a = -G\K b + 2i 1 A a 5 1 = % -A a A a , 

5G\ = G a c G c b + c ace c bde A c A d +A a A b + ^A c A c 5 a b . (2.37) 

Note that the variation of G a b indeed leaves invariant the cubic norm N[Z] for an anti- 
commuting A a . 

We now consider the equations of motion for the scalar fields parametrizing the sym- 
metric space These are expressed in terms of the corresponding Maurer-Cartan form 

v~ 1 dv = P + B, (2.38) 

which decomposes accordingly into the coset component P defining the scalar momenta, 
and the SI2 ® 04 component B defining the pulled back spin connection. In components, 
the scalar momenta are defined as 

w = -dU - ^e 2U *dcu, S a = -efdt* , Z = e u Z{kF), Z a = e u ' Z a {-kF) , (2.39) 

where we introduce some shorthand notations that will be used for the remainder of the 
section. Analogously, we give the components of B along s[2 

p(B) = -\e 2U *duj - ±Q, X(B) = e u Z{*F) , (2.40) 

its components along 04 

l{B) = -le 2U *doj + ±Q, A a (B) = e u Z a (*F) , (2.41) 

and finally G a b (B) defines the 64 valued traceless component of the pulled back spin con- 
nection 011M4: 

G\(B) = efdjel dp - e^e^df - , (2.42) 

where Q is the pulled back Kahler connection ([2,8]) . 2 

These formulae given, one can straightforwardly compute the equations of motions of 
the scalar and vector fields respectively as coming from the equation of motion and Bianchi 
identity on P, as follows 

d B *P = 0, d B P = 0, (2.43) 

where d B stands for the covariant derivative on the coset. For instance, the components of 
the equation of motion for P are 

d B *w = -d-kdU- ^e 4U duj*du+ e 2U Re[Z(F)Z(*F) + Z a (F)Z a (*F)] = 0, 
d B *£ a = -V-kefde + e 2U (2Z(F)Z a ( y *F) + c abc Z b {F)Z c {*F)) = , (2.44) 



To prove that G a b{B) is indeed traceless, one can use fl2.20| ) to show that c c bcd = ^^^j* and 
substitute (EyJ) in c abc d^c abc = 0. 
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and 



d B * Z = e u [DZ{F) - eld? A Z a (F) - ie 2U * dto A Z[F)\ = , (2.45) 
d B * Z a = e u (DZ a {F) - e a% dP A Z(F) - c a6c efdf A Z C (F) + ie 2U *du> A Z a (F)) = , 



where (2.27) and the standard special geometry identities ( 2.10| ) were used. In the following 



we will not analyse these equations directly, but will rather employ arguments based on 
the nilpotency of P for extremal solutions to obtain equivalent first order equations that 
can be solved directly. 

3. Nilpotent orbits and first order systems 



In this section we generalise the formalism developed in [45] to arbitrary nilpotent orbits 
of G3, with a specific emphasis to the ones describing non-BPS black holes. The basic 
observation is that the only regular stationary solutions of J\f = 2 supergravity with a flat 
three-dimensional base metric are such that the momentum P is nilpotent. This implies 
in particular that P can be written in a basis of generators e a which lie in a nilpotent 
subalgebra of 33. Such a nilpotent subalgebra is always associated to a semi-simple element 
3 h of sfe © g4 such that 

he a := [h,e Q ] = p a e a , 1 < p a < n , (3.1) 

where n defines the maximal possible eigenvalue of adh in 33. This implies for instance the 
equation 

n 

l[(h-i)P = 0, (3.2) 

i=i 

which defines a first order constraint between the components w, Z, Z a and S a of P. In 
order to be consistent with the equations of motion and the Bianchi identity ( [2.43 ), the 
covariant derivative of the generator h must satisfy 

n n i—l n n i—l 

n o*-j)d B hH(h-k)AP=o, n (h-*hn(h-^p=o. (3.3) 

i=l j=i+l k=l i=l j=i+l k=l 

These equations are satisfied if d B h also lies in the nilpotent algebra defined by h, or 
equivalently that 

n 

]J(ad h - i)d B h = . (3.4) 

i=i 

In general, one can always choose the generators h such that only its components A and 
A a do not vanish. Equation ( |3.4| ) can be viewed as first order equations for these auxiliary 
components, which can be solved to determine their evolution in space in terms of the 
physical fields. As mentioned above, equation fl3.2|) defines first order equations for the 
physical fields which contain these auxiliary components A and A a and determine dU + 



3 Semi-simple means that it is in the conjugation class of an element of the Cartan subalgebra. 
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^e 2U *du) and e^dt 1 in terms of e u Z(*F) and e u Z a (*F), plus some possible constraints on 
the latter if the dimension of the coset component of the nilpotent algebra defined by h is 
strictly less than 2n v + 2. 

For BPS solutions one has A a = and A = e ia , where the phase a defines the covari- 
antly constant spinors as in ||, and ( |3,4j) is equivalent to the equation 

da + Q- ^e 2U *dui = 0. (3.5) 

The non-BPS solutions with vanishing central charge at the horizons are described in a 
similar fashion, with A = and a normalised rank one A a (i.e. c abc A b A c = and A a A a = 1) 



1 45]. We will now discuss the specific examples of the nilpotent orbits associated to the 
systems describing respectively composite and single centre non-BPS black holes with a 
non- vanishing central charge at the horizon. 

Composite nilpotent elements 

The composite non-BPS solutions admit a scalar momentum P which lies in the posi- 
tive grade components of the graded decomposition of the coset component 2 <S> M 2riv+2 , 
associated to an element h of $34 that leads to the decomposition 

04 = (R n «) ( - 2) e (six e g 5 ) (0) e (R n -) (2) , (3.6) 

for Q4 itself and 

2 ® R 2 ""+ 2 2 ( " 3) (2 ® R n ") (_1) (2 ® R n ") (1) © 2 (3) , (3.7) 

for the coset component, i.e. P € (2 ® R n ") (1) © 2 (3) . Such an element h can always be 
chosen to be Hermitian (i.e. to lie in 04 (u(l) ©64)) so that it is realised for A (h) = fi a , 
where f2 a satisfies 

Af[n]n a = ^ c abc n b n c . n a n a = 3 . (3.8) 

Equivalently, f2 a is in the U(l) x K± orbit of the Jordan algebra identity. 
More explicitly, one finds the following action on the coset component 

hw = n a z a hz = n a z a hz a = n a w + c abc n b ^ c hs a = n a z + c abc n b z c . (3.9) 

Considering the grade three part of P, from the equation [h, P (3) ] = 3P (3) one obtains the 
solution 

Z (3) = A/p]w (3) Z a(3) = f] a w (3) S a(3) = N[n]ti a w^ , (3.10) 

for an arbitrary w (3) . Similarly, from the equation [h, = for the grade one part, 
one obtains the solution 



w 



(1) - n a z a(1) z {1) = -M[np a z^ s a(1) = c abc n b z^ -N[n}n a n b z ( b 1) , (3.11) 



for an arbitrary Za 



(1) 
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Considering a general linear combination of these two solutions one concludes that Z 
and Z a are arbitrary, whereas w and S a are determined as 

w = ^(n a z a - Af[n]z) , n a = c abc n b z c + ^ti a (z-N[n]n b z b ), (3.12) 

which are the explicit first order relations for the scalar momenta. 

These contain the auxiliary components jV[f2] and O a of h in (|3.8[) , which can be viewed 
as defining a very small vector R of unit mass (i.e. I±(R) = dI±(R) = d 2 I^(R)\ Slj A = and 
\Z(R)\ = 1) through 

Z(R) = N[Q] Z a (R) = Q a . (3.13) 
The flow equations for these fields are given by ( |3.4| ), which in this system reduces to 

[h,d B h] = 2d B h . (3.14) 

Using the explicit form of B ( 2.4C )- (|2.42j ) and the first order constraint fl3.12|) , one computes 
the components of d B h as 

7 (d B h) = hm[n a z a ], 



A a (d B h) = Z a (dR) + N[tt]e al dt J + c abc Sl l 'efdtf - -e 2U * dun o 



z a {dR) + Re[Af[n]n b eldt l }z a (R) - ( z a + M[n]c abc n b z c - n a n b z 



G a b (d B h) = c ace c bde [n d z c - n c z d ) + n b z a - ti a z b - ^im[n c z c ] , (3.15) 

where we explicitly separated the terms depending on the derivative of the vector R. It 
is now straightforward (though cumbersome) to compare ( prill ) with the ab ove relations, 
using that 

c ace c bde n d (Af[n]c cfg n f z 9 - n c n f z f ) - n a (Af[n]c bcd n c z d - n b n c z c ) 
= -c ace c bde n c z d + n b z a , (3.16) 

which follows by ( 2.20| ). The result is that ( |3.14j ) is satisfied provided that 

dR = -Re[AT[n]n 6 eJdt*]fl , (3.17) 

which implies that there exist a constant symplectic vector R such that 

„ R /x 

R = s— . (3.18 

\Z(R)\ 

We conclude that the generator h is in this case determined by a constant very small 
projective vector R and the scalar fields such that 

0.19) 

One can now return to ( |3.12j ), which becomes a first order flow equation for the scalars e u , 
•kduj and dt l in terms of the gauge fields and the constant vector R. 
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Single centre nilpotent elements 

The composite non-BPS system above can be consistently reduced to a system describing 
the single centre class of solutions. The associated graded decomposition consists in break- 
ing furthermore the s[2 algebra by introducing a non-compact generator h* defined such 
that its only nonvanishing component is A = e ia . The action of this generator follows from 
( |2~33|) as 

h*w = e" ia Z, KZ = e ia w, KZ a = e -ia E a , h*£ a = e iQ Z a , (3.20) 

and, by its own, it would define the BPS system. The single-centre non-BPS solution is 
defined in the positive grade component of the graded decomposition associated to the 
generator ^(h + h*) 

2 x R 2n "+ 2 r ( ~ 2) e (leR"")'- 1 ' e (R n * ©rj (0) e (k®e"") (1) ©i (2) , (3.21) 

where h is the generator that defines the composite system above. One can straightfor- 
wardly compute that the solution ( 3.1C| ) decomposes into [h*, P]_ 3) ] = ±P± } according to 

wf = ±e- ia M[Q]wf (3.22) 

and these two solutions define the grade 1 and the grade 2 sing lets in (|3.21| ). On the other 
hand, the grade 1 component of (|3.21| ) in R nv corresponds to the solution ( |3.11| ) satisfying 
moreover IK,P^ } } = Pj 15 , or explicitly 

Z a + (1) = e~ ia {c abc VL b Z£l - N[n\ti a ti b Z£l) . (3.23) 

Summing up the two solutions, one obtain that the single centre non-BPS momenta satisfy 
( 3.12 ) for Z and Z a constrained to satisfy the phase dependent equation 

z a - M[n\n a z = e~ ia (c abc n b z c + n a (z- M[n\n b z b )) , (3.24) 

which represents a constraint on the physical degrees of freedom that is necessary to reduce 
to single centre solutions. 

This is expected for a single centre solution, since the element fi a is defined by its 
overall phase and the angle K4/K5 ( being the maximal compact subgroup of G5, and 
therefore the stabilizer of £l a in -K5), that is n v real parameters in total. The constraint 
( 3.24 ) defines precisely n v real equations, such that one can think of it as determining fl a 
in terms of the central charge Z, its derivatives and the 'BPS phase' a. 

On the other hand, when viewed as a constraint on the charge vector it is simple to see 
that ( p. 24 ) reduces its components by half. Defining the combination K a = Z a — J\f[0]O, a Z, 
the constraint becomes 

K a = e~ ia {c abc tt b K c - Af[n] Cl a Cl b K b ) = i(K a ) . (3.25) 

In the right-hand-side we defined the operation l, which is an anti-involution 

i(i{R a )) = e- ia (c abc n b i(K c ) -M[fi]n a n b i(K b )) 
= c abc n bCcde n d K e - n a n b R b 

= K a , (3.26) 
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where we used (|3.8| ) in the last equation. In what follows we will elaborate on these points 
of view of the constraint ( p. 24 ), in connection to the various aspects of the solutions. 

We close this section by giving the analog of the consistency condition (|3.14| ) for the 
generator h*. In this case, only the 5(2 components are important, and using fl2.40| ) one 
computes that 



p(d B K) = -2lm[e- ia Z] 



\(d B h* 



De ia + | e 2U * dw e ia 



(3.27) 



where the covariant derivative D is the Kahler covariant derivative in (2.9), consistent with 
the unit Kahler weight of the phase e ia . Imposing that this is a grade two element of the 
5(2 algebra, i.e. that 

2p{d B K) = i(e ia \(d B K) - e- ia X(d B K)) , 2\{d B K) = 2ie ia p{d B K) , (3.28) 
turns out to be equivalent to the purely imaginary condition 

e -ia jj^ia j + i e 2U + duJ = _ 2 ■ j m [g-fa^r] ^ (3 29) 

which fixes the phase a in terms of the physical degrees of freedom through 

da + Q + ±e 2U *dw = -2Im[e- ia Z}. (3.30) 

Since we also have the phase JV[fi], it is natural to define the Kahler invariant phase 
e ia N[tt]. From JOUD we find that 

darg[A/"[fi]] -Q + Im[A/"[^]O a efdf] = , (3.31) 

and using moreover (3.12) one obtains that 

d(a + arg[A/p]]) + ^e 2U * du 

+ lm[Af[ti}n a e?dt i + e- ia Af[n](dU - ie 2U *duj - Af^nae^df)] =0. (3.32) 

Through ( 3.12[ ), ( |3.24| ) also implies a reality constraint on the scalar field momenta 

Mp)e1dt + e ia Np] (-c abc n b e cT df l + AA[ft]M b e te dP) 

= 9, a (dU - ie 2U *daj + e ia Af[n}(dU + §e 2C/ *dw)) (3.33) 

and in particular 



//[npaefdt* + e ia Q, a e al dt~ l = 3{dU - \e 2V * doo + e ia M[Q](dU + \e 2V * doo 
Therefore one obtains finally 

d(a + arg[A/p]]) = e 2U * duo - 2Im[e ia Af[ti} (dU + |e 2C/ * dui)] 
\ - cos(a + arg[AT[0]])) f e 2U * duj 



(3.34) 



tan i 



q+arg[Af[n]] 



dU , (3.35) 
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which can be integrated to 



_ e -2U 



*du = d — m = dM . (3.36) 

We conclude that the angular momentum is given in terms of a harmonic function M 
dual to the Kaluza-Klein vector, which is a known characteristic feature of single centre 
solutions [|29[j. The phase e' ia M[0,] is determined in this way in terms of spacetime fields as 

(-M + ie~ 2U ) 2 
= ( J "J • (3-37) 



e 



+ M 2 



4. Duality covariant form of the non-BPS black hole solution 

In this section we solve explicitly the flow equations in the first order system describing 
non-interacting non-BPS black holes discussed in the last section. We will see that in 
this case the vector fields w and uj carrying the electromagnetic charges and the angular 
momentum are simply sourced by harmonic functions, although the vector fields satisfy a 
quadratic constraint such that they only depend on n v + 1 harmonic functions, instead of 
2n v + 2 in the BPS system ||. Incidentally we exhibit that this first order system reduces 
to a linear system of differential equations. This is a necessary step towards the explicit 
solution of the non-BPS composite system describing interacting non-BPS black holes. 
After presenting the procedure of integrating ( [3.12 ) combined with the reality constraint 



( 3.24 ), we briefly discuss the physical properties of the solutions. 



4.1 Integrating the first order equations 



The starting point is the solution of the nilpotency condition (|3.12|) , written explicitly as 



a first order system for the scalars and the metric degrees of freedom 

dU + l -e 2U *duj = --e u (n a Z a {-kF) -Af[Q]Z(*F)) , 

-efdt* = c abc Q b e u Z C (*F) + ±ti a e u (Z(*F) - N[n]tt b Z b (*F)) , (4.1) 

where F is the spatial component of the field strengths defined in Q2.26P . For later reference, 
we give the inverse relations for the field strengths 

e u Z(*F) = ^A/p] (dU - l -e 2U * dwj - -^df , 

e u Z a (*F) = -c abc n b e^df + ^MppAdf - ~n a (dU - % -e 2U * doo) , (4.2) 
where we used the short-hand notation f2j = e^O, a . 
Electromagnetic scalar potentials 

In order to solve this system, we combine the information on the derivative of R given by 
(3.17) with (4.2) to construct the gauge field momenta as 



d( = 2 e 2U Re [Z(*F)V + Z i (*F)D i V] 

= d (e u Re [Af[ti]V - WDiV]) + \e u Im [Mppidt] R + \e 3U *duR, (4.3) 
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where we made extensive use of the special geometry identities ( 2.10| ). The first term is 
manifestly a total derivative, whereas the others are along the very small vector R, and 
must therefore combine into the derivative of a single function. This implies the existence 
of a function M such that 

e u Im [Af[Q] n t dt l ] + ±e 3C/ * doj = M e 3U Re [jV^fijdi*] - d(M e 3U ) . (4.4) 

And indeed, in the 'single centre' system one shows using ( ET37D in <$M) that M is the 
function that determines the angular momentum in ( 3.36| ), It follows that the gauge field 
momenta take the form 

C = e^Re [Mp[V - WDiV] - ^e 3U MR , (4.5) 
with the corresponding central charges given by 

Z(C) = \ e u (1 + ie 2U M)M[Q] , Z a {Q = % - e u (1 + * e 2U M)Q a , (4.6) 
for later reference. 

Using ( [4.1| ) and the structure of the ( |4.5| ), one shows that one vector is always trivial, 
because 

(R,Q = -2e u \Z(R)\, (4.7) 
whereas, taking the imaginary part of (fO|) one gets 



e 2U * duo = — (R, *F) , (4. 

2\Z(R)\ 



and using ( 2.26j ), one gets therefore that 

(R, dw) = . (4.9) 

Note that this property does not require the reality constraint ( 3.24| ) and is also valid for 
composite non-BPS solutions [44]. 

The linear system 

One can now combine (^1]) and (^1]) to disentangle the term proportional to -kdu in the 
definition of the scalars, such that ( f4.lQ becomes 



+\e w {-M + ie- 2U )M[np a * du , (4.10) 



-e\dt = c abc n b e u Z c (*dw) + \SFe v (Z(*dw) - N[n]n b Z b (*dw)) 

l 

2 



which, using the constraint ( 3.24 ), can be rewritten as 



e»dt* = e u e ia (Z a (*dw) - N[n]n a Z(*dw)) -M[n]n a (dU- -e 2U *doj) . (4.11) 
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Applying the same procedure on fl4.2|) , one obtains the the inverse relations of ( 4,11| ) 
for the central charges Z{dw) and Z a (dw). The charge vectors dw can then be straightfor- 
wardly constructed with the result 



dw = 21m [-Z{dw)V + Z a (dw)D a V] 

= -2 e- u lm [e~ ia (- * dU + \e w duj)V + e'^dfDiV + \ fiTl] 



(4.12) 



where we used the shorthands 



// = -^e 4U d(e~ w + M 2 ) + e~ ia N[£l] (*dU - ^e 2U duj - Afp]^ * df) . (4.13) 

At this stage we have exhausted the constraints implied by the existence of the constant 
R, and it is important to find a constant vector Darboux conjugate to R in order to be 
able to decompose conveniently dw. This is indeed possible, using equations (3.17), ( |3.33j ), 
(|3736b , Q) and Q, which allow one to show that 



R* = \Z(R)\- 1 Re Y 3 Af[Q] V + \Y\ 2 YSt D { V 
Y = (1 + i e 2U M) , 

is constant, in the following way 

1 



(4.14) 
(4.15) 



dR* = -2R* Re[Af[n]nidf] + ^Im[3Y 2 AA[0]V + Y{2Y - Y)Q i D i V]d(e 2U M) 

\ z i R )\ 

Re[Y 2 ti T df l V + Y(YAf[n]df + Y ^c^VL^dt^) AV] 



\Z(R)\ 



2R* R&lMp^idt] + ^lm[Y 2 M[n]V + |y| 2 ff D t V] d{e 2U M) 

\Z(R)\ 



\Z(R)\ 



Re[(Y 2 Af[n]V + |y| 2 ff AV) (AA[0]adf)] 



0. 



(4.16) 



This vector is indeed mutually non-local with R, since their inner product, (R*,R) 



and is also very small, i.e. its central charges satisfy (3^8), as a consequence of its definition. 
One can now project dw along this new vector to find 



{R*,dw) 



*d(e- u \Z(R)\- 1 \Y\ 2 ) = -*dV, 



(4.17) 



where we defined the distinguished harmonic function V, whose pole will carry the linear 
combination (R ,T) of the physical charges. Note that V can equivalently be defined as 



V 



-U 3 



'4\Z{R 



*M2 



\Z{R)\ 



(4.18) 



and because the central charge of a very small vector is nowhere vanishing in moduli 
space, it follows that for a regular extremal solution, the function V is strictly positive and 
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(R*,T) < 0. Combining this with fl3"l7|) , (gg ) and ( pl|) we can determine the combination 
Sljcft* in terms of V, the metric components and the phase e ia as 



Re [Af[n] n 4 dt* 



\Y\~ 2 d\Y\ 2 



e- ia £lidt l = -FRe[AT[0]Oi(ff] 



(4.19) 



We can now use this information to further simplify the expression for dw by adding 
multiples of the two constant vectors with appropriate coefficients. It turns out that the 
most suggestive form is obtained by subtracting a multiple of R and arranging that the 
sign of the i?-component is flipped. This clearly reflects the situation one encounters in 
the known explicit solution defined through the almost-BPS system [26, and indeed 
one finds that 



dw 



(R*,dw) 
(R*,R) 



R 



-21m [e- u e~ %a ((- * dU + \e 2U dw)V + dt B&)\ 

+2Ve u Jm (ie~ ia V) *d(f ) + *d(f ) R* 
= -2Im*D(e- u e- ia V) + *d(f) R* , 

where we defined the modified covariant derivative 

D(e- u e- ia V) =[d + i{Q + da + \e 2U * dw - e 2U d(f ) V)] {e~ u e~ ia V) . 



(4.20) 



(4.21) 



The form of this covariant derivative is exactly such that the corresponding composite 
connection is trivial by use of ( |3.30| ), ( |4.2| ) and ( |4. 19 ) 



Q + da + 

It then follows that ( [4.20 ) takes the form 



3 2U d(f )V 



0. 



21m {e~ u e- ia V) - (f ) R* 



dw-2"-?^R = -*d 

(R*,R) 

which implies that the vector fields are defined in terms of harmonic functions as 

dw = -kdH , 



(4.22) 



(4.23) 



(4.24) 



such that for instance (R*,H) 
are given by 

2Im(e- u e-" x 



-V . Using this back in (4.23) one finds that the scalars 



V) 



M 



(R*,H) 



R* . 



(4.25) 



Note that this result only depends on harmonic functions, where the angular momentum 
harmonic function only appears through its ratio with V, as in [21], and the harmonic 
functions H control the electric and magnetic charges. However these harmonic functions 
are not all independent as in the BPS system, but are subject to algebraic constraints that 
reduce them to n v + 1 independent functions, as we show in the next paragraph using the 
constraint ( |3.24 ). This constraint is rather non-linear because of its dependence on the 
scalar fields, and we will now rewrite it as a quadratic constraint in the harmonic functions 
themselves. 
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Quadratic constraint 

For this purpose will make use of the quartic invariant of TV = 2 supergravity coupled to 
a symmetric scalar manifold, which reads [^] 

h(T) = ^t MNp QT M T N T P T Q 

= -(P^l) 2 +±qvN[p\ - ip°M[q] + c ijkI Pp k c ilm qi q m , (4.26) 

in terms of the cubic norm. Here, t MNp Q is a completely symmetric tensor, and M, N, . . . 
are symplectic indices that encompass both the upper and the lower components in (|2.2j). 
The absolute value of this expression is known to determine the entropy of static black 
holes for any value of the charges. We also define the lift of symplectic indices through the 
symplectic form 

r*-(£), r- = (7), (4,7) 

such that 

(r l5 r 2 ) = r 1M rf = -r 2M rf . (4.28) 



We know already from ( |4.9[) that the harmonic functions dH = *dw must be symplectic 
normal to the vector R, and because the integration constant of % along R can always 
be reabsorbed in a redefinition of the integration constant of the function M in ( 4.25| ), we 
can assume without loss of generality that 

(H,R) = 0. (4.29) 

To rewrite the constraint (13.241) we will consider the vector 



h M CH,R*) ^ Mum N A * N = \ tMNPQ KNU P R* Q , (4.30) 



where R* is the small vector defined in ( 4.14 ) as the magnetic dual to R. To compute the 
decomposition of this vector in terms of its components linear in T-L M itself and the small 
vectors R M and R* M , it will be useful to observe the following consequence of symplectic 
invariance 



§ {q , P ) = - iZ (%, - f) , §- a = iZ a {%, - f) , (4.31) 

which will allow us to compute ( (4.30| ) in the complex basis starting from the alternative 
expression for the quartic invariant of a vector in terms of its central charges 

h = {ZZ - Z a Z a f - c eab Z a Z h c ecd Z c Z d + 4ZA/"[Z] + 4ZA/"[Z] . (4.32) 

Note that, while the scalars appear explicitly in all terms, the complete expression can be 
shown to be scalar independent and equal to ( [4.26 ). Using the above properties and the 
reality constraint ( |3.24 ) on dM. = -kdw, one can compute that 

t MNPQ QvUp = 2 q^ u) g u)H M _ 4 <i?*,W(i?*,W fl M (4 33) 

(R* , R) 
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The interested reader can find an outline of the derivation in appendix [B[ This constraint 
can be integrated to the same constraint on the harmonic functions themselves, up to 
possible constants which do not a priori need to satisfy the constraint. Nevertheless, using 
J4.25 ), and substituting the expressions ( 3.37D , ( 4.14 ), ( 4.15 ) and ( [4.17 ) one computes that 



Z(R - \VR) = e- u+ia + yZ(R*} 



-e- u lMM + l -Me u ^MM 7 



Z a (U-\VR) = yZ a (R* 



--Me u Ytt a . 



(4.34) 



It follows that R — \VR satisfies the constraint ( 3.24 ), which implies that R does as well. 
Since R must satisfy (3.24), it follows that the integration constants in R also satisfy ( 4.33| ) 
and one obtains that 



1 



I± M (n,R*) = {R*,R)R 



M 



{R*,R) 



(4.35) 



We can now use the fact that both R and its derivative satisfy the above constraint, to 
show that they must necessarily lie in a Lagrangian subspace. Indeed, for any two vectors 
1^2 satisfying the constraint (3.24), one can show that their inner product ( p. 15 ) in the 
complex basis can be written as 



<ri,r 2 



2Im 



-z(ri)z(r 2 ) + k ^a(ri) n b z b (r 2 ) 



(4.36) 



which vanishes upon requiring that Ti^ are mutually local with R, as expressed by (|B.5| ). 
It then follows that 

(H,dH) = 0. (4.37) 

We conclude that the poles of R must be mutually local charges, and such solution cannot 
describe interacting black holes. Indeed, in Appendix [A] we show an example of this prop- 
erty in a specific duality frame to obtain the the n v + 1 harmonic functions parametrising 
the relevant Lagrangian subspace. 



Scaling factor and moduli 

This concludes our analysis. For any very small vector R one can construct explicit solu- 



tions, after first using the asymptotic moduli in order to determine R* from ( 4.14 ), which 
by definition will be constant. One must then solve the algebraic equation ( |4.35[ ), which 
determines the allowed harmonic functions R. The scalars and the metric scale factor can 
be obtained by solving (4.25) in the standard way [53|. For instance 



-4(7 



h(R 



->VR 



fR*) 



(4.38) 



Using the property that R* is very small, it follows that the terms of order three and four 
in M vanish and the term of order two simplifies according to j54fl 



U MNPQ R* M R* N T P r Q = -(R*,T)< 



(4.39) 
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Considering moreover that ( 4.35|) is satisfied for T~L — \ VR and contracting this constraint 
with U-\VR, one obtains that the term of order one in R* vanishes as well. Therefore 



-4C7 



h(H - \VR) - *C {R*,H - \VR) 2 = h(H - \V R) - M 1 



(4.40) 



Using again ( [4.34 ) one computes that the component H. — \ VR which mutually commutes 
with R* satisfies 



T Y 2 M[n] , 



-u. 



2\Y\ 



z a (n-\vR) = ^ 

which permits to compute that 

h(n - \vk) = o . 



Y 2 )n a 



(4.41) 



(4.42) 



It follows that I/£H — hV R) is linear in V, and since the factor — h specifically switches 



the sign of the R component one concludes that 

e~ w = -I 4 (H)-M 2 . 
One obtains in the same way that the moduli can be expressed as 



1 dl 4 
2&Hi 



(n) + %e-* u VR* + (M + ie 



1 dh 
"2 



§^CH) + ie- 2U VR° + (M + ie~ 2U ) [% 



o 



v 



-R*°) 



(4.43) 



(4.44) 



after several simplifications. The metric is then given by ( |2.24j ) with oj as in (|3.36 ), whereas 



the gauge fields are given by (2.26) with <^ as in ( |4.5[ ) and dw given by (4.24) above. 
4.2 Physical properties 

In the preceding subsection, we have treated the constant vector R as defining the system, 
and the constraint ( |3.24j ) or ( |4.35 ) as a restriction on the physical charges for a given R. 
Physically, it is however more natural to define a solution from the asymptotic moduli 
f^, the electromagnetic charge T and angular momentum J. Considering the asymptotic 
central charge ^(r)^ and its Kahler derivative ZaiT)^ (which we will refer to as the 
'central charges' for simplicity), one can indeed define the asymptotic VL aoo as the unique 
solution of ( 3.24 ) for which Oqo is determined such that there is no NUT charge, i.e. 

Im[O a Z a (r) -MmZiT)]^ = , (4.45) 



as it is done in [44]. Indeed, it follows from ( |4.1| ) that ( 4.45 ) defines the NUT charge and 
that 



1 



n a z a (r)-M[Q]z(r))\ oo = M A 



(4.46) 



is the non-BPS fake superpotential at spatial infinity, i.e. the ADM mass. If one does not 
fix the gauge for the U(l) x K4 gauge invariance, the 'central charges' do not depend on 
the flat directions in moduli space, and it follows that f2 aoo then does not depend on the 
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flat directions either. However, the constant vectors R and R which are defined from 
^aoo upon action of the asymptotic moduli through ( |3.19| ) and (4.14) do. 4 Therefore the 



asymptotic 'central charges' altogether with the constant vector R contain the information 
about the flat directions. 

To understand this property, let us discuss the stabilizers of R and T in G4. It is known 
that the stabilizer of a very small vector as R is [55| 



G 5 x R™ 11 C G 4 , (4.47) 

whereas the stabilizer of the electromagnetic charges with a strictly negative quartic in- 
variant ^(r) < is G5 C G4. 5 However, as we exhibit in Appendix [A], only the compact 
subgroup C G5 of the stabilizer of the charges leaves the very small vector invariant. 
It follows that the action of the non-compact generators which generate the flat directions 

G 5 /K 5 c G 4 /^4 , (4.48) 

act faithfully on R. Moreover, one can show that the condition that both R and R* are 
very small, altogether with equation ( [4.35| ) 

-l' M (F,R*) = (R*,T} T M -2 {R¥,T)2 R M , (4.49) 

entirely determines these small vectors up to an overall rescaling in terms of the electro- 
magnetic charges and n v — 1 parameters parametrizing the flat directions. We prove this 
in Appendix [A] in a specific duality frame. 

Considering a single centre solution carrying charges T and angular momentum J 

r cos 9 

U = h + - , M = m + J^- , (4.50 

the scalar fields on the horizon take the form 



• jfftCr) + ^s^f- + (j cos e + is) (p + l s {R , 



|^(r) + 2*S^ + (Jcos^ + iS) (TO + z5-g^ 



(4.51) 



where 



S = V~h(T) - J 2 cos 2 9 . (4.52) 

As is clear from ( |4,51| ) , the attractor values of the moduli are not entirely determined by the 
electromagnetic charges T and the angular momentum J, but depend on the asymptotic flat 
directions through the small vectors R and R , in general. This formula ( 4.51| ) generalises 
the rotating attractor formula derived in [28, ^] in specific duality frames. 



4 Note that Y<x> is determined by the asymptotic 'central charges', because Ooo is. 

5 This can easily be checked for a D6 very small vector which only non- vanishing charge is p°, because the 
latter is clearly left invariant by the five-dimensional duality group G5 and the n v T-dualities. In the same 
way, a D0-D6 charge with only non- vanishing qo and p°, is clearly left invariant by G5 only, and admits a 
negative quartic invariant. 
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Although the scalar fields are not entirely determined by the electromagnetic charges 
and the angular momentum, it follows from ( 4, 43] ) that the horizon area only depends on 
the electromagnetic charges and the angular momentum as expected (2^, |24], 57] 



A = 4ny/-U(T) - J 2 . 
The same formula implies that the ADM mass is determined as 



Ma 



where we assumed that asymptotically 



e~ 4U \ 

c I oo 



1 dhQX) M 

4 dh M 



- J 4 (h) - m 2 = 1 . 



(4.53) 



(4.54) 



(4.55) 



Although this formula may suggest that the ADM mass is linear in the charges, one must 
note that its explicit expression in terms the asymptotic 'central charges' is generally a non- 
rational function of the latter [fl7f| , due to the fact that the integrating constants h are not 
entirely parametrized by the asymptotic moduli alone. The situation is similar, although 
simpler, for BPS black holes, for which the asymptotic moduli are entirely determined in 
terms of the integrating constants h of the harmonic functions dual to the electromagnetic 
vectors, but the reverse is not true, as the constants h are parametrized by the asymptotic 
moduli and the phase of the asymptotic central charge. Indeed, the ADM mass of a BPS 
black hole is not a linear function of the central charge 



1 dh(h) M 
4 dh M 



\Z(T)\ C 



(4.56) 



due to the presence of precisely this phase. For the non-BPS solutions, the constants h are 
also not entirely determined by the asymptotic moduli, since they depend explicitly on the 
asymptotic 'central charges' through the phase of the central charge as well as the small 
vectors R and R* . Moreover, the asymptotic moduli are not entirely determined in terms 
of the constants h either in this case, as is clear from the asymptotic value of ( [4.44 ) 



2 dhi ' 



(4.57) 



This expression includes explicitly the small vectors that depend on the fiat directions, 
and therefore do not affect the mass formula, but includes also the constant m which 
parametrizes the phase of e _JQ A/'[0]|oo, and therefore depends explicitly on the asymptotic 
'central charges'. 

The regularity of the solution requires that e~ 2U is everywhere strictly positive, and 
the absence of closed time-like curves outside the horizon moreover requires that e 



-4U 



> 



' J sin 9 



2T^) . The latter condition implies the former, and it reads 



-h(H) > m + 



\J\ 



(4.58) 
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This condition is clearly satisfied at spatial infinity because of ( [4.55 ), and at the horizon 
this requires the usual regularity condition 



-7 4 (r) > J 2 , (4.59) 
which is necessary for the horizon area ( f4.53j ) to be well defined. 

5. Conclusion 

In this paper we have given a detailed exposition of the first order systems underlying 
the composite non-BPS system of multi-centre black holes in N = 2 supergravity in four 
dimensions with a symmetric very special Kahler geometry. Upon imposing a reality con- 
straint on the system of equations, we restricted to the single centre class, which includes 
all extremal under rotating solutions with one centre and multi-centre generalisations with 
mutually local charges. Making use of this constraint we were able to explicitly integrate 
the flow equations for the single centre class for the vector multiplet scalars in a manifestly 
duality covariant way. Here, we discuss some of the implications of our results. 

The solution we obtain for the single centre class, being manifestly duality covariant, 
allows for general moduli at infinity and arbitrary charge vectors, without the need of 
dualising a specific seed solution. We stress the presence of an additional (constant) very 
small vector R and its magnetic dual R* in the solution for the moduli, in addition to the 
standard vector of harmonic functions describing the charges. This vector arises in the 
definition of the flow equations for the full composite non-BPS class and therefore plays a 
central role in our considerations. This is quite different from the squaring of the action 
in the standard fake superpotential approach for single centre solutions [H], 15, 19, 17, 



18, 58], which is based on a function of scalars and physical charges only. In view of the 



fact that our explicit solution ( [4.25 ) allows to construct a function driving the flow that 



contains R, R* along with the charges and scalars, it remains an interesting open problem 
to understand the relation between the two formulations. 

From a physical point of view both R and R* are integration constants for the scalar 
equations of motion once the charges are fixed. Indeed, for a single centre solution with 
given charges it is known that not all scalars take part in the flow from infinity to the 
near horizon region, but particular combinations are frozen to arbitrary constant values 



throughout spacetime, the so called flat directions |27|, |59j, |60|, |55j, [yj. We have shown that 
the ambiguity in defining R from the electromagnetic charges is precisely parametrized by 
the flat directions in moduli space, as expected in order to describe single centre solutions 



explicitly. Addition of more centres with charges such that the constraint ( 4.35 ) is satisfied 
lifts the flat directions, since a unique R is fixed in terms of the charges in the generic case. 

The possible microscopic interpretation of the general single centre under-rotating 
solution remains unclear at the moment. A number of approaches have been proposed for 



the microscopic construction of extremal non-supersymmetric black holes, see e.g. 61, p2| 



63]. From this point of view, the flat directions appear as geometric moduli or background 



fluxes that can take arbitrary values, see e.g. |l3|, 33] 
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A natural future direction is the construction of the generic solution in the composite 
non-BPS class. Since the flow equations are again characterised by the very small vector R, 
it is clear that some of the structures found here will remain relevant in the more general 
case. The almost BPS class should also admit a similar description. 

Finally, the recent results in J64| , 65] indicate that a non-extremal deformation or a lift 
of our flow equations to five and/or six dimensional supergravity would be very interesting 
to explore, in connection to the over rotating branch. 



Acknowledgement 

We thank Alessio Marrani for useful discussions. This work was supported by the French 
ANR contract 05-BLAN-NT09-573739, the ERC Advanced Grant no. 226371, the ITN 
programme PITN-GA-2009-237920 and the IFCPAR programme 4104-2. 



A. Under rotating seed solution 

In this Appendix we present the known rotating seed solution in a specific duality frame 
1 29], as a convenient pivot to draw intuition for the general solution. In this case the 
electromagnetic vector fields satisfy 

*dw = ——=dV , -kdw 1 = —=dU , dwi = dw° = . (A.l) 
v 2 V 2 



The scalar fields take the form 

,- —M + ie' 2U 



L\ (A.2) 



N[L\ 
and the metric 

e~ 4U = VN[L] - M 2 , *du = dM . (A.3) 
In this duality frame the constant small vectors are 

R= (o, 0; 2V2, o) T , R* = (V2,0;0,o) T , (A.4) 



which satisfy indeed 

Z(R*) = iY 3 Z(R), D l Z(R*) = -i\Y\ 2 YD l Z(R), (A.5) 

where we used the definition fl4.15| ) and M is the harmonic function in (A.3). One can 
straightforwardly solve ( |4.35| ) in terms of these vectors using that 



1 d 2 i 



4 o + Qop° = -qop° - QiP 1 



2 dqodp' 

1 ^ 2 ^4 i ijk 

1 d 2 h 

q qi = -2q qi 



2 dpidp 



1 d 2 L 



1 q q + 2q 2 = . (A.6) 



2 dp dp 
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Because %q = —-j^V by definition, the third line implies that %i = 0, and then T-L° = 0, 
which was already implied by (4.29). Therefore we find the consistent solution 



(A.7) 



It is then straightforward to check that e iU is indeed equal to ( 4.43| ) and that the moduli 
are equal to Q4.44D 

. (M + ie~ 2U )^L i 



1 M[L] 



-V- 



9^(M + %e~™)y/2 

-M + ie~ 2U ri 
L 



M + ie~ 2U M[L] 
Let us now consider the general solution of ( [4.49D for given charges 

r= (o,^;^, o) T , 



(Ai 



(A.9) 



associated to such a solution. For this purpose, we will parametrize the two small vectors 
in terms of n v + 1 parameters as 

R = c (-Af[f\, -hc ijk fjfk ; 1, ~f) T , R* = c* (l, e l ■ -AT[e], \c i]k eie k ) T , (A.10) 

which can in general be infinite, as long as R and R* are themselves finite in the limit. Let 
us solve ( (4.49|) as an equation for these two small vectors. The component along Rq of this 
equation implies that 



fc\2 



(<?o - lc ijk p l e J e k ) 



and substituting this into the Ri component, one obtains 

fi-- 



(A.ll) 



M[e] 3 k 

2 C-ijkV P 



2g 2 



(A.12) 



Substituting these equations and the R l component of ( 4.49| ) inside its R° component one 
finally obtains 

(A.13) 



\c ijk p i p j e k = —Af[e]Af\p] 



1 

qo 



which altogether imply that ( 4.49 ) is satisfied identically. 

One observes that the small vectors ( A.10| ) can then be obtained from the ones in ( |A.4[ ) 
by a duality transformation that acts on the scalar fields as 

t\e) =f + e l + ^ k c ]lv c kqr e l pPtH r - - t Vc jkl e k p l + 0(e 2 ) , (A.14) 

at first order (if one fixes the irrelevant constants c = 2\/2 and c* = V2). One straight- 
forwardly computes that such transformations leave the charges T invariant. In order for 
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the moduli to be well defined, the stabilizer of the charges p l in G5 must be its maximal 
compact subgroup K§. Therefore it follows that the n v — 1 e*'s which parametrize the 
solution of ( 4.49| ), also parametrize the moduli space of flat directions G5/K5 C G4/K4, 



and alternatively, that equation ( [4.35 ) uniquely determines the very small vectors in terms 



of the electromagnetic charges, up to n v — 1 parameters associated to the flat directions. 
B. Duality invariant constraint 

In this appendix we give an outline of the derivation of ( 4.35| ) from the reality condition 



( 3,24 ) for a vector J that is mutually local with the vector R, i.e. (J, R) = 0. To this end, 



we compute the vector in ( 4.30| ) 



l' 4 M ( J, J,R*) = l -t MNp Qj N J P R* Q , (B.l) 

where R* the small vector symplectic dual to R defined in ( |4.14| ). Note however that, since 
the computation is homogeneous with respect to all vectors, we will rather use 

Z(R*) = e^M^ [Cl] Z a [R*] = e^M^ [n]Q a , (B.2) 



here, or in other words we replace R* by its associated vector of mass one. We now proceed 
to compute the components of the derivative ( |4.30|) in the complex basis, as in( 4.31 ), 



starting from the expression ( 4.32 ). Writing Z, Z a for Z(J), Z a (J), one obtains that 



- e-^M^\n\c ah si h c cde z d z e - 2e'fAr 1 i[n]c abc z b c cde n d z e 

- e ! fM^[ri}Q a (ZZ - Z b Z b ) -2Z a Re[2e^N^[fl} - N^[£l]tt a Z a ] . (B.3) 
Using the reality constraint (|3.24 ) one eliminates all the terms in c abc Cl b c cde Z d Z e ^ and using 



it again on the resulting expression one eliminates the terms in c abc Z b c cde VLdZ e . The final 
expression is then 

1 01 3 

--J- (J, J 7 R*) = -2iZ a lm\Ze^ 1 M^[n]-e^M^[n]n a Z a } - Sl a (ef N^p]ZZ 

2 oZ a 1 \ 

+(e^Afa [fl]Z - e^N 1 * [n]n a Z a + e'^Af^ p]Z) (Af[£l]Z - n a Z a + e ia Z + 2e" ia Z)) . 

(B.4) 

In order to simplify the second term we need to use the property that (J, R) = 0. This 
permits to solve for 

_ 3e^ A/j [Q] Z + e%Af* [O] (Af[ti]Z + 2Af[Q]Z) 
e"3\Af5[fi] + e--M2[Q] 

which, when used in ( |B.4| ) leads to 



1 dl 



2dZ a 



(J,J,R*) = -2iZ a Im[Ze'^M^[Q] - et M^[n]Q a Z a ] 



2 ilm[e ajVzffi]] 



+ 4n a lm[e-* a Zy ic^'J - (B.6) 
Re[e-M2[fl}} 2 
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In the same way, one computes the Z component of the derivative 



-^(J,J,R*) = 2iZIm\Ze^Afhn] - e*f N^[£l]£l a Z a ] 
2 oZ 



- 4M\n]lm[e- ia Z} 

Refe^AAf^]] 2 

In order to interpret these two equations, let us compute that 

(R*,J) = 2Im[e-^AA2[0]Z -e- l -?N^[n}ti a Z a ] 

.Im[ e -f.A^[Q]] 2 



-4Im[e~ tQ Z} 



Re[e^A/-2[f2]]2 



and 



(R*,R) = 2Jm[e- S ^Af^[n]Af[n] - e _ ^.A/i[n]n a n o ] 
= -8Im[e-^M^[n}} 3 . 

Using these identities we conclude that (B.6) and flB.71 ) combine into 



1 dh 



2 8J M 



(J,J,R*) = {R*,J)J 



M 



'(R*,R) 



(B.7) 



R 



M 



(Bi 



(B.9) 



(B.10) 



Using the homogeneity of this equation in R and R* , one can write it for the constant 
vectors R and R* such that this equation defines a quadratic algebraic equation in J, as 



claimed in section 4.2. 



References 

[1] S. Perrara, R. Kallosh, and A. Strominger, J\f = 2 extremal black holes, Phys. Rev. D52 



(1995) 5412-5416, |iep-th/9508072 |. 
[2] A. Strominger, Macroscopic entropy of 'N ' = 2 extremal black holes, Phys. Lett. B383 (1996) 



39-43, ]hep-th/9602iii l. 

[3] S. Ferrara and R. Kallosh, Supersymmetry and attractors, Phys. Rev. D54 (1996) 1514-1524, 



[ hep-th/9602136 ]. 

[4] K. Behrndt, D. Lust, and W. A. Sabra, Stationary solutions of M — 2 supergravity, Nucl. 



Phys. B510 (1998) 264-288, frep-th/9705169 |. 
[5] F. Denef, Supergravity flows and D-brane stability, JEEP 08 (2000) 050, |hep-th/0005049(| 



[6] G. Lopes Cardoso, B. de Wit, J. Kappeli, and T. Mohaupt, Stationary BPS solutions in 
Af = 2 supergravity with R 2 interactions, JEEP 12 (2000) 019, ]hep-th/0009234(|. 



[7] J. P. Gauntlett, J. B. Gutowski, C. M. Hull, S. Pakis, and H. S. Reall, All supersymmetric 
solutions of minimal supergravity in five dimensions, Class. Quant. Grav. 20 (2003) 



4587-4634, |iep-th/0209114 



[8] J. P. Gauntlett and J. B. Gutowski, General concentric black rings, Phys. Rev. D71 (2005) 



045002, [hep-th/0408122l 



-27- 



A. Castro, J. L. Davis, P. Kraus, and F. Larsen, String theory effects on five-dimensional 
black hole Physics, Int. J. Mod. Phys. A23 (2008) 613-691, 



arXiv: 0801. 1863 



A. Ceresole and G. Dall'Agata, Flow equations for non-BPS extremal black holes, JHEP 03 



(2007) 110, [hep-th/0702088 



G. Lopes Cardoso, A. Ceresole, G. Dall'Agata, J. M. Oberreuter, and J. Perz, First-order 
flow equations for extremal black holes in very special geometry, JHEP 10 (2007) 063, 



arXiv: 0706. 3373 . 



A. Ceresole, G. Dall'Agata, S. Ferrara, and A. Yeranyan, First order flows for Af = 2 extremal 



black holes and duality invariants, Nucl. Phys. B824 (2010) 239-253, farXiv : 0908 . lliofl 
E. G. Gimon, F. Larsen, and J. Simon, Black holes in super gravity: the non-BPS branch, 



JHEP 01 (2008) 040, |arXiv:0710.4967[ ]. 

D. Gaiotto, W. W. Li, and M. Padi, N on- super symmetric attractor flow in symmetric spaces, 



JHEP 12 (2007) 093, [arXiv:0710. 16381 



L. Andrianopoli, R. D'Auria, E. Orazi, and M. Trigiante, First order description of black 
holes in moduli space, JHEP 11 (2007) 032, |arXiv:0706.0712| . 



L. Andrianopoli, R. D'Auria, E. Orazi, and M. Trigiante, First order description of D = 4 
static black holes and the Hamilton- J acobi equation, Nucl. Phys. B833 (2010) 1-16, 



[ arXiv: 0905. 3938[ | . 

G. Bossard, Y. Michel, and B. Pioline, Extremal black holes, nilpotent orbits and the true fake 
superpotential, JHEP 01 (2010) 038, 



arXiv: 0908. 1742 



A. Ceresole, G. Dall'Agata, S. Ferrara, and A. Yeranyan, Universality of the superpotential 
for d = 4 extremal black holes, Nucl.Phys. B832 (2010) 358, |arXiv:0910.2697| . 



J. Perz, P. Smyth, T. Van Riet, and B. Vercnocke, First- order flow equations for extremal 



and non-extremal black holes, JHEP 03 (2009) 150, |arXiv:0810. 1528J , 

S.-S. Kim, J. Lindman Hornhmd, J. Palmkvist, and A. Virmani, Extremal solutions of the S 3 
model and nilpotent orbits o/G 2 (2)> JHEP 1008 (2010) 072 



arXiv: 1004.5242 



P. Galli, K. Goldstein, S. Katmadas, and J. Perz, First- order flows and stabilisation 
equations for non-BPS extremal black holes, JHEP 1106 (2011) 070, ]arXiv: 1012.4020[ . 



D. Rasheed, The rotating dyonic black holes of Kaluza-Klein theory, Nucl. Phys. B454 



(1995) 379-401, |hep-th/9505038 l. 

T. Matos and C. Mora, Stationary dilatons with arbitrary electromagnetic field, Class. Quant. 



Grav. 14 (1997) 2331-2340, |hep-th/9610013 l. 

F. Larsen, Rotating Kaluza-Klein black holes, Nucl. Phys. B575 (2000) 211-230, 



| hep-th/9909102 l. 

T. Ortin, Extremality versus supersymmetry in stringy black holes, Phys. Lett. B422 (1998) 



93-100, flhep-tri/9612142 |, 

R. R. Khuri and T. Orti'n, A Nonsupersymmetric dyonic extreme Reissner- Nordstrom black 



hole, Phys. Lett. B373 (1996) 56-60, [hep-th/9512178] 



P. K. Tripathy and S. P. Trivedi, N on- super symmetric attractors in string theory, JHEP 03 



(2006) 022, [hep-th/0511117l 



- 28 - 



D. Astefanesei, K. Goldstein, R. P. Jena, A. Sen, and S. P. Trivedi, Rotating attractors, 



JHEP 10 (2006) 058, fhep-th/0606244fl . 

I. Bena, G. Dall'Agata, S. Giusto, C. Ruef, and N. P. Warner, Non-BPS black rings and black 



holes in Taub-NUT, JHEP 06 (2009) 015, [nrXiv: 0902 .4526 



G. Dall'Agata, S. Giusto, and C. Ruef, U-duality and non-BPS solutions, JHEP 1102 (2011) 
074, |arXiv:1012.4803|| . 



K. Goldstein and S. Katmadas, Almost BPS black holes, JHEP 05 (2009) 058, 



arXiv: 0812. 4183 



G. Bossard and C. Ruef, Interacting non-BPS black holes, Gen.Rel.Grav. 44 (2012) 21-66, 



arXiv: 1106.5806 



I. Bena, S. Giusto, C. Ruef, and N. P. Warner, Multi-center non-BPS black holes - the 
solution, JHEP 11 (2009) 032, |arXiv:0908.2121 |. 

I. Bena, S. Giusto, C. Ruef, and N. P. Warner, A (running) bolt for new reasons, JHEP 11 



(2009) 089, ]arXiv: 0909. 2559[ | . 

I. Bena, S. Giusto, C. Ruef, and N. P. Warner, Supergravity solutions from floating branes, 



JHEP 03 (2010) 047, [ |arXiv:0910. 186C| ]. 

N. Bobev and C. Ruef, The Nuts and Bolts of Einstein-Maxwell Solutions, JHEP 1001 



(2010) 124, | |arXiv: 0912. 0010| | . 

P. Breitenlohner, D. Maison, and G. W. Gibbons, Four- dimensional black holes from 
Kaluza-Klein theories, Commun. Math. Phys. 120 (1988) 295. 

M. Giinaydin, A. Neitzke, B. Pioline, and A. Waldron, BPS black holes, quantum attractor 



flows and automorphic forms, Phys. Rev. D73 (2006) 084019, [ |iep-th/0512296| ]. 

E. Bergshoeff, W. Chemissany, A. Ploegh, M. Trigiante, and T. Van Riet, Generating geodesic 



flows and supergravity solutions, Nucl.Phys. B812 (2009) 343-401, [ |arXiv: 0806. 2310 1. 
G. Bossard, H. Nicolai, and K. S. Stelle, Universal BPS structure of stationary supergravity 



solutions, JHEP 07 (2009) 003, JarXiv: 0902. 4438 1. 

G. Bossard and H. Nicolai, Multi-black holes from nilpotent Lie algebra orbits, Gen. Rel. 



Grav. 42 (2010) 509-537, |arXiv : 0906 . 1987| ]. 

P. Pre, A. S. Sorin, and M. Trigiante, Integrability of supergravity black holes and new tensor 



classifiers of regular and nilpotent orbits, JHEP 1204 (2012) 015, JarXiv: 1103.0848( |. 

W. Chemissany, P. Giaccone, D. Ruggeri, and M. Trigiante, Black hole solutions to the 
F^-model and their orbits (I), arXiv : 1203 . 633g| . 

G. Bossard, Octonionic black holes, 



arXiv: 1203.053C. 



G. Bossard, 1/8 BPS black hole composites, 



arXiv: 1001.3157 



R. Kallosh, N. Sivanandam, and M. Soroush, Exact attractive non-BPS STU black holes, 
Phys. Rev. D74 (2006) 065008, thep-tri/0606263| . 

B. de Wit and A. Van Proeyen, Potentials and symmetries of general gauged M = 2 
supergravity: Yang-Mills models, Nucl. Phys. B245 (1984) 89. 

B. de Wit, P. G. Lauwers, and A. Van Proeyen, Lagrangians of M = 2 supergravity - matter 
systems, Nucl. Phys. B255 (1985) 569. 



- 29 - 



A. Ceresole, R. D'Auria, and S. Ferrara, The symplectic structure of M = 2 supergravity and 



its central extension, Nucl. Phys. Proc. Suppl. 46 (1996) 67-74, ||hep-th/9509160 l. 

M. Giinaydin, G. Sierra, and P. Townsend, The geometry o/A/"= 2 Maxwell-Einstein 
supergravity and Jordan algebras, Nucl. Phys. B242 (1984) 244. 

S. Ferrara and S. Sabharwal, Quaternionic manifolds for type II superstring vacua of 
Calabi-Yau spaces, Nucl.Phys. B332 (1990) 317. 

S. Ferrara, E. G. Gimon, and R. Kallosh, Magic supergravities, M = 8 and black hole 



composites, Phys.Rev. D74 (2006) 125018, [hep-th/0606211 



B. Bates and F. Denef, Exact solutions for super symmetric stationary black hole composites, 



JEEP 1111 (2011) 127, |hep-th/0304094 



S. Ferrara, A. Marrani, E. Orazi, R. Stora, and A. Yeranyan, Two-center black holes 
duality-invariants for STU model and its lower-rank descendants, 



arXiv: 1011.5864 



S. Ferrara and A. Marrani, On the moduli space of non-BPS attractors for J\f = 2 symmetric 



manifolds, Phys. Lett. B652 (2007) 111-117, |arXiv : 0706 . 1667| ], 

S. Ferrara and R. Kallosh, On JV = 8 attractors, Phys.Rev. D73 (2006) 125005, 



| hep-th/0603247 l. 

A. Sen, Black hole entropy function and the attractor mechanism in higher derivative gravity, 



JHEP 09 (2005) 038, |hep-th/0506177|l 



A. Ceresole and S. Ferrara, Black holes and attractors in supergravity, arXiv: 1009.4175 



S. Ferrara and M. Giinaydin, Orbits of exceptional groups, duality and BPS states in string 



theory, Int. J. Mod. Phys. A13 (1998) 2075-2088, |iep-th/9708025 |. 

S. Bellucci, S. Ferrara, M. Giinaydin, and A. Marrani, Charge orbits of symmetric special 



geometries and attractors, Int. J. Mod. Phys. A21 (2006) 5043-5098, |hep-tti/0606209|1 
R. Emparan and G. T. Horowitz, Microstates of a neutral black hole in M theory, 



Phys. Rev. Lett. 97 (2006) 141601, [tiep-th/0607023| 



A. Dabholkar, A. Sen, and S. P. Trivedi, Black hole microstates and attractor without 



super symmetry, JHEP 0701 (2007) 096, [tiep-th/0611143 



E. G. Gimon, F. Larsen, and J. Simon, Constituent model of extremal non-BPS black holes, 



JHEP 07 (2009) 052, [ |arXiv:0903.0719| ]. 

P. Galli, T. Ortin, J. Perz, and C. S. Shahbazi, Non-extremal black holes of N — 2, d = 4 



supergravity, JHEP 1107 (2011) 041, |arXiv: 1105.3311 1. 

I. Bena, M. Guica, and W. Song, Un-twisting the NHEK with spectral flows, 



arXiv: 1203.4227 



- 30 - 



